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ABSTRACT 

Cosmic metal enrichment is one of the key physical processes regulating galaxy forma¬ 
tion and the evolution of the intergalactic medium (IGM). However, determining the 
metal content of the most distant galaxies has proven so far almost impossible; also, 
absorption line experiments at 2 : 6 become increasingly difficult because of instru¬ 

mental limitations and the paucity of background quasars. With the advent of ALMA, 
far-infrared emission lines provide a novel tool to study early metal enrichment. Among 
these, the [C ii] line at 157.74 ^m is the most luminous line emitted by the interstellar 
medium of galaxies. It can also resonant scatter CMB photons inducing characteristic 
intensity fluctuations (A///cmb) near the peak of the CMB spectrum, thus allowing 
to probe the low-density IGM. We compute both [C ii] galaxy emission and metal- 
induced CMB fluctuations at z ^ 6 by using Adaptive Mesh Refinement cosmological 
hydrodynamical simulations and produce mock observations to be directly compared 
with ALMA BAND6 data (r'obs ~ 272 GHz). The [C ii] line flux is correlated with 
Muv as log(Fpeak/MJy) = —27.205 — 2.253 Muv — 0.038 M^v Such relation is in very 
good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 
2015) of Muv < —20 galaxies. We predict that a Muv = —19 (Muv = —18) galaxy 
can be detected at 4cr in ~ 40 (2000) hours, respectively. CMB resonant scattering can 
produce ~ ±0.1 /iJy/beam emission/absorptions features that are very challenging to 
be detected with current facilities. The best strategy to detect these signals consists 
in the stacking of deep ALMA observations pointing fields with known Muv — —19 
galaxies. This would allow to simultaneously detect both [C ii] emission from galactic 
reionization sources and CMB fluctuations produced by 2 ~ 6 metals. 

Key words: cosmology: cosmic microwave background - infrared: general - galaxies: 
high-redshift - galaxies: intergalactic medium 


1 INTRODUCTION 

In ACDM cosmologies^, structure formation is a hierarchi¬ 
cal bottom-up process (i.e. Press & Schechter 1974). Dwarf 
galaxies are expected to be abundant and to represent the 
first (2 ± 10) efficient metal factories of the Universe (e.g. 
Madau et al. 2001; Ferrara 2008). These galaxies effectively 
pollute their surrounding intergalactic medium (IGM), since 
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^ In this work we assume a ACDM cosmology with total matter, 
vacuum and baryonic densities in units of the critical density 
Oa = 0.727, = 0.228, Ulf, = 0.045, Hubble constant Hq = 

100 h km s“l Mpc“^ with h = 0.704, spectral index n = 0.967, 
(Tg = 0.811 (Larson et al. 2011). 


- compared to larger galaxies - dwarfs have a shallower po¬ 
tential well and because their smaller size allows multiple 
SN to coherently drive the outflows (e.g. Ferrara & Tolstoy 
2000 ). 

In spite of the impressive progresses produced by deep 
optical/IR surveys (Dunlop 2013; Madau & Dickinson 2014; 
Bouwens et al. 2014) in identifying galaxies well within the 
Epoch of Reionization, very little is known about the metal- 
licity and other properties of these systems, including feed¬ 
back (e.g. Dayal et al. 2014) and interactions with their en¬ 
vironment (e.g. Barnes et al. 2014). 

The situation is slightly better for what concerns the 
IGM, where metal enrichment is typically studied by mea¬ 
suring the abundance of heavy elements at different cosmic 
times with quasar (QSO) absorption line spectroscopy (e.g. 
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Songaila 2005; Ryan-Weber et al. 2009; Becker et al. 2009; 
Simcoe et al. 2011; D’Odorico et al. 2013) and - more re¬ 
cently - with gamma-ray bursts soft X-ray absorption (e.g. 
Campana et al. 2010; Behar et al. 2011; Campana et al. 
2015). These observations show that the diffuse IGM is en¬ 
riched at metallicity Z ^ 10~^'®Zq at any overdensity (A) 
and redshift {z) probed so far (Meiksin 2009). As the gas left 
over by the Big Bang is (virtually) metal-free, understand¬ 
ing how and when the first metals were produced holds the 
key of many structure formation processes. 

As a consequence of the paucity of standard luminous 
lighthouses prior to reionization (Barkana & Loeb 2001; Cia- 
rdi & Ferrara 2005), the enrichment at 2 6 cannot be effi¬ 

ciently studied with absorption line experiments. Addition¬ 
ally, at high -2 the detection of most metal tracers becomes 
increasingly difficult. For example, C iv metal absorption 
line detections at high -2 are rare since the C iv doublet shifts 
into the near-infrared. In this spectral regions observations 
become challenging because of increased OH emission from 
the sky and severe telluric absorption. 

The above situation is going to change very soon 
thanks to the advent of the Atacama Large Millimeter- 
submillimeter Array (ALMA). This observatory gives us ac¬ 
cess to far-infrared (FIR) and sub-mm spectroscopy, where 
metal lines carrying information about the energetics of the 
interstellar medium and outflows are found. FIR emission 
lines are a powerful tool to study high -2 galaxies, since - un¬ 
like Lya - they are affected neither by the increasing IGM 
neutral fraction at 2 6 nor by the presence of dust. 

In particular, the [G ii] (^^ 3/2 ^^ 1 / 2 ) line at 

157.74 /rm is the brightest FIR emission line, i.e. it can ac¬ 
count for ~ 1% of the total infrared luminosity of galaxies 
(e.g. Grawford et al. 1985; Madden et al. 1997). Before the 
advent of ALMA, [G ii] observations at 2 4 were limited 

to QSO host galaxies and to rare galaxies with extreme star 
formation (~ IO^Mq yr“^)(e.g. Maiolino et al. 2005; Garilli 
& Walter 2013; Gallerani et al. 2012; Cicone et al. 2015). The 
unprecedented resolution and sensitivity of ALMA offer the 
unique opportunity to search for [G ii] emission from nor¬ 
mal star forming galaxies (~ 1 — IOMq yr“^) at high -2 (e.g. 
Maiolino et al. 2015; Capak et al. 2015). These observations 
yield key information about the stellar age, gas metallicity, 
star formation, and also on the environments in which these 
early galaxies formed (e.g. Gonzalez-Lopez et al. 2014; De 
Looze et al. 2014). Recently, the [G ii] emission arising from 
high -2 star forming galaxies has been modelled through nu¬ 
merical simulations (e.g. Nagamine et al. 2006; Vallini et al. 
2013; Tomassetti et al. 2015; Vallini et al. 2015) and ana¬ 
lytical models (e.g. Gong et al. 2012; Munoz & Furlanetto 
2014). 

In addition to tracing the ISM of galaxies, FIR reso¬ 
nant lines correspond to transitions that are excited via a 
resonant scattering process of Gomic Microwave Background 
(CMB) photons on heavy element atoms present in the IGM 
(e.g. de Bernardis et al. 1993; Maoli et al. 1996). These res¬ 
onant transitions are available for essentially all the most 
abundant species, e.g. C, N, O, Si, S and Fe, in the mid-IR 
and FIR wavelength ranges (e.g. Basu et al. 2004). 

As a result, IGM metals can produce GMB spectral dis¬ 
tortions and spatial fluctuations that can be used to extract 
unique information on the cosmic metal enrichment process, 
like, for example, spatial distribution maps of a given metal 


species. This experiment is conceptually similar to QSO ab¬ 
sorption line studies, but it uses the GMB as a background 
source: the method has the enormous advantage that every 
pixel in a sky map can act as a source against which metal 
lines may appear either in emission or absorption. 

Various works (e.g. Maoli et al. 2005; Basu 2007; 
Schleicher et al. 2008) have proposed GMB fluctua¬ 
tions/distortions as a tool to trace the smooth distribution 
of metals and molecules in the post-recombination Universe 
(2 3> 10). These fluctuations typically affect the first GMB 
multipoles - at I ^ 10^, large angular scales {d ~ 1 deg) 
- and they will be hopefully seen by Planck successors 
(e.g. Hernandez-Monteagudo et al. 2006; Ghluba 2014). On 
smaller scales - 0 ~ 10” {I Si 10®) - CMB fluctuations in the 
FIR can be used to map the enriched IGM surrounding the 
first galaxies (2 10). 

Here we present a model based on detailed Adaptive 
Mesh Refinement hydrodynamical cosmological simulations, 
that allows us to simultaneously compute both the [C ii] line 
emission from early galaxies and the expected level and sta¬ 
tistical properties of CMB fluctuations arising from inter- 
galactic C II resonant scattering. 

The paper is organized as follows. In Sec. 2, we describe 
our model, introducing the adopted high -2 cosmic metal en¬ 
richment simulation (Sec. 2.1). We postprocess the simu¬ 
lation to account for [C ii] galaxy emission (Sec. 2.2) and 
CMB fluctuations (Sec. 2.3), and analyze the properties of 
the simulated CMB fluctuations in Sec. 3. Then, we con¬ 
struct and analyze mock ALMA observations of [C ii] galaxy 
emission and CMB fluctuations (Sec. 4). Finally, we propose 
and discuss an observational strategy for the detection of 
[C ii] galaxy emission and metal-induced CMB fluctuations 
(Sec. 5). Conclusions are given in Sec. 6. 


2 MODEL 

2.1 Cosmological simulations 

We start from the suite of cosmological hydrodynamical 
metal enrichment simulations presented by Pallottini et al. 
(2014a, P14 hereafter). Using a customized version of the 
Adaptive Mesh Refinement code ramses (Teyssier 2002), 
P14 follows the evolution of a (10 Mpc h~^)^ comoving cos¬ 
mic volume from 2 = 199 until 2 = 4 with 512® dark matter 
particles. The gas evolution is tracked on a corresponding 
number of coarse grid cells, with 4 additional levels of re¬ 
finement based on a Lagrangian mass threshold-based crite¬ 
rion. This set-up allows us to reach a maximum resolution 
of 1.22 h~^ kpc for the gas in the densest regions. 

Star formation and supernova feedback are included via 
customized subgrid prescriptions; P14 accounts for stellar 
yields based on population synthesized models (Salvador! 
et al. 2008). The simulated galaxy sample reproduces the 
observed cosmic star formation rate (Bouwens et al. 2012; 
Zheng et al. 2012) and stellar mass densities (Gonzalez et al. 
2011) evolution in the full redshift range 4 ^ 2 10. 

In P14, we track the progressive metal enrichment of the 
gas, initially of primordial (BBN) composition, and its ther¬ 
mal history by accounting for heating and cooling processes 
(Theuns et al. 1998). The gas ionization state is regulated 
by an external, redshift-dependent ionizing UV background 
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Figure 1. Left panel: UV luminosity function, >I>uv at z = 6. Red circles represent the simulated 'I’uv; blue squares denote the 
observed "huVi inferred from a composite collection of HST datasets (Bouwens et al. 2014); black dashed lines are the Schechter fit 
(Bouwens et al. 2014) to the observations. Right panel: [C ll] flux as a function of Mjjv at z = 6 (blue small circles); the same relation 
binned in AMpv = 1 intervals is shown by red circles. The red errorbars correspond to the r.m.s. variance within the magnitude intervals; 
The red dashed line is the best fit to the relation. Also plotted are the data from recent high-z observations from Maiolino et al. (2015); 
Capak et al. (2015); Willott et al. (2015) with yellow squares, green pentagons and violet hexagons, respectively. 


(Haardt & Madau 1996, 2012, UVB) produced both by stars 
and QSOs. At z ~ 6 , the observationally allowed rage for 
the UV photoionization rate is ~ (0.9 —3.6) x 10“^^s“^ (e.g. 
Wyithe & Bolton 2011). We have verified that the gas ther¬ 
modynamical evolution is only marginally affected by UVB 
intensity variations within this range. 

To compute the ionization state of the various atomic 
species, we postprocess the simulation outputs using the 
photoionization code CLOUDY (Ferland et al. 1998). Simi¬ 
lar to P14, we account both for the UVB intensity at 912 
A (J,/), and the density (n), kinetic temperature (T*,), and 
metallicity {Z) of the gas. 

We adopt the following classification for the baryonic 
gas: (i) the interstellar medium (ISM) defines highly over- 
dense gas (n/rimean = A > 10®); (ii) the true intergalactic 
medium (IGM) is characterized by A 51 10; (iii) the circum- 
galactic medium (CGM) represents the interface between 
the IGM and ISM (10 < A < 10®). Pallottini et al. (2014b) 
have shown that the adopted classification is consistent with 
the one proposed by Shull (2014), based on the distance r 
between the gas and the center of the galaxy. Specifically (i) 
the ISM is located within the galaxy virial radius (rvir), (ii) 
the IGM at r/rvu ~ 5 and (iii) the CGM at 1 r/rvir ~ 5. 

In P14, we showed that galaxies develop a stellar mass- 
metallicity (M*-Z*) relation by z = 6. For M* > IO^Mq, 
the stellar metallicity [Z^) increases with increasing stellar 
mass, while for M* IO^Mq, Zi, is constant. This means 
that, while massive galaxies are able to retain most of their 
metals, low mass galaxies are prone to metal ejection be¬ 
cause of their swallower potential well. This indicates that 
at high-z low mass galaxies are the main driver of metal 
enrichment in the IGM (e.g. Ferrara 2008). In agreement 
with previous numerical studies (e.g. Dayal et al. 2013), the 
relation shows little evolution from z = 6 to z = 4, 
for galaxies with M* ~ 10^. Our results are consistent with 
the metallicity evolution suggested by recent observations of 
3 < 2 < 5 star forming galaxies (Troncoso et al. 2014). 


We find that (i) most of the gas resides in the IGM, 

(ii) the denser CGM contains about 15% of the baryons and 

(iii) the ISM accounts only for a small fraction of the to¬ 

tal mass (~ 7%). On the other hand, at any given redshift 
the enriched gas is mostly found near the metal production 
site, i.e. in the ISM; less than 10% in mass of the produced 
metals can reach the IGM/CGM. In particular, at z = 4, a 
A-Z relation is in place: the IGM shows an uniform distri¬ 
bution around Z ~ 1O“® ®Z 0 , in the CGM Z steeply rises 
with A up to 1 O“®Z 0 and the ISM has Z ~ 1 O“®Z 0 . A 
considerable fraction ( 50%) of the enriched IGM/CGM 

is in an hot state {Tk ^ 10‘^ ®K), meaning that most of the 
carbon is in C iv. Our analysis shows that G iv absorption 
line experiments can only probe ~ 2 % of the total produced 
carbon. 

P14 results are in agreement with previous numerical 
studies in terms of the baryon thermodynamical state (e.g. 
Rasera & Teyssier 2006; Cen & Chisari 2011), the evolution 
of the metal filling factor (e.g. Oppenheimer et al. 2009; 
Johnson et al. 2013) and the A-Z relation (e.g. Gnedin & 
Ostriker 1997; Oppenheimer et al. 2012). Additionally, the 
preliminary analysis of synthetic spectra extracted from the 
simulation is consistent with recent observations of metal 
absorption lines (D’Odorico et al. 2013). 

As a further consistency check, we test the UV lumi¬ 
nosity function (4>uv) extracted from the P14 simulation 
against observations of z = 6 galaxies. For each galaxy in our 
simulation, Muv is calculated with STARBURSt99 (Leitherer 
et al. 2010 ), using Z* and the stellar age as input parame¬ 
ters. In the left panel of Fig. 1, we plot the simulated and 
observed "Fuv with red circles and blue squares, respectively. 
The observations are taken from a composite collection of 
HST datasets (Bouwens et al. 2014), and we additionally 
plot the Schechter fit to the observations with a black dashed 
line. The simulated 4?uv well matches the observed UV lu¬ 
minosity function for —20 Muv ~ — 16. Additionally, the 
simulated 'Fuv is in fair agreement with the Schechter fit 
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extrapolated at —16 ^ Muv ~ — 12. Note that, for the sim¬ 
ulated $uv, the bright end is around Muv ~ ~19) where 
the large scatter is due to the limited simulated volume. 


2.2 [C ii] emission from the ISM of high-z galaxies 

Starting from the P14 simulations, we compute the 
[C ii] emission at 157.74 ^m due to the ^i- 3/2 —> '^Pi /2 for¬ 
bidden transition of ionized carbon in the ISM of galaxies. 
The [C ii] line is the dominant coolant of the galaxy in¬ 
terstellar medium, and it can be collisionally excited under 
conditions present in different ISM phases, e.g. in the cold 
and warm neutral medium (CNM, WNM), in high density 
photodissociation regions (PDRs), and even in the ionized 
gas (Tielens & Hollenbach 1985; Wolfire et al. 1995, 2003; 
Abel 2006; Vallini et al. 2013). 

In our cosmological hydrodynamic simulations, we can¬ 
not resolve sub-kpc scales, that are typical for different 
ISM structures (CNM, WNM and PDRs). To overcome this 
problem, we account for [C ii] emission^ from the ISM of 
high- 2 ; galaxies by adopting the model presented in Vallini 
et al. (2013) and updated in Vallini et al. (2015) (V15 here¬ 
after). V15 predicts the [C 11 ] emission from CNM, WNM 
and PDRs within individual galaxies as a function of their 
metallicity (.Z*) and star formation rate (SFR). The V15 
model is based on the radiative transfer simulation of a high- 
2 galaxy; it includes a subgrid treatment for computing the 
thermodynamical equilibrium of the diffuse neutral medium 
(on ~ 60pc scales); it allows to localize molecular clouds 
within the galaxy (on ~ 1 — 10 pc scales) and to compute 
the corresponding [C ii] emission by means of the PDR code 
UCL_PDR (Bell et al. 2005). 

According to the V15 model^, the [C ii] galaxy lumi¬ 
nosity (I/cii) depends on the galaxy SFR and metallicity 
through the following fitting formula: 

log (Lcii) = 7.0 -b 1.2 log(SFR) -b 0.021 log(Z*) (1) 
-b 0.012 log(SFR) log(Z*) - 0.741og^(^*), 

where Ten, SFR and Z^, are in units of Lq, Mq yr“^ and Z©, 
respectively. The V15 model has been used to interpret high- 
2 : LAEs and LBGs [C ii] observations (e.g. Gonzalez-Lopez 
et al. 2014; Ota et al. 2014), and V15 provides predictions 
consistent with [C ir] observations of local metal-poor dwarf 
galaxies (De Looze et al. 2014). We compare the V15 pre¬ 
dictions with the work by Munoz &; Furlanetto (2014) (see 
their eq. 8 ): for a SFR = lO^M 0 /yr and a solar metallicity, 
according to the V15 model Leu — 2.5 x lO^L©, that is a 
factor 2 ::^ 5 higher that the one predicted by Munoz &; Furlan¬ 
etto (2014), i.e. Lcii ^ 5 x 10®L©. Given the uncertainties 

^ In the present work, we focus on [C ii] emission from 2 ~ 6 
galaxies. However, we note that the V15 model can be readily 
extended to compute the emission of fine structure lines from 
other heavy elements, such as, e.g. [N ii] at 122 pm and [O i] at 
63 pm (Vallini et al. 2013). 

^ The V15 model assumes that H ii regions provide only a neg¬ 
ligible contribution to the [C ii] emission (Gonzalez-Lopez et al. 
2014, in particular see Fig. 8). However, the relative fraction of 
[C ii] emission arising from H ii regions is not clearly known. As 
an indicative value, Vasta et al. (2010) observe a median contri¬ 
bution of ~ 30%. For a more detailed discussion, we address the 
interested reader to Yue et al. (2015). 
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Figure 2. [C ll] CMB spectral distortions/fluctuations from a 
simulated line of sight (l.o.s.) at 2 = 6. Left panel: intensity 
fluctuations (AA/A) as a function of frequency (u, lower axis) 
and velocity {v, upper axis). Absorption and emission features are 
highlighted in red and blue, respectively, and the corresponding 
peaks are labelled with Ai and Ei. Note that integration in the 
shown bandwidth {Au ~ 0.1 GHz) yields A/// ~ —7 X 10“^. 
Right panel: overdensity (A, top), metallicity {Z, center) and 
peculiar velocity {vp, bottom) for the [C ll] metal patch respon¬ 
sible for the fluctuations (|A7i.//i.| > 0) as a function of co¬ 
moving distance (kpch~^). In the overdensity plot our A based 
IGM/CGM definition is marked with a dashed black line. In the 
peculiar velocity plot tip = 0 is marked with a dashed black line. 
The origins of both the velocity (left panel) and distance (right 
panel) axis are shifted for displaying purpose. 


on the local [C ii]-SFR relation, 6 x 10® Tcii/T© ~ 10^®, 
we find that both models are consistent with observations, 
though the V15 predictions provide a better match. To com¬ 
pute Lcii for galaxies in the P14 simulation we use eq. 1 , 
where .Z* and SFR are the mass-averaged values derived 
from the simulation (see also Yue et al. 2015). 

We calculate the relation between [C ii] and UV emis¬ 
sion. In the right panel of Fig. 1, we show the [C 11 ] flux as 
a function of Muv for the simulated galaxies at 2 = 6 and 
for various 2 = 5 — 7 ALMA observations (Maiolino et al. 
2015; Capak et al. 2015; Willott et al. 2015). The simulated 
relation is shown as a scatter plot, binned in AMuv = 1 in¬ 
tervals. The plot shows a clear correlation between [C ii] flux 
and UV magnitude for the simulated galaxies, with a Spear¬ 
man correlation coefficient of ~ —0.78. The best fit formula 
for the relation is 

log(F’peak/MJy) = -27.205 - 2.253 Muv - 0.038 Myv ( 2 ) 

Using eq. 2 we can slightly extrapolate the simulated rela¬ 
tion to Muv < — 20 , and directly compare it with recent 
ALMA observations (right panel of Fig. 1), i.e. a [C ii] de¬ 
tection in a 2 ~ 7 LBG (Maiolino et al. 2015, yellow square), 
observations of several “normal” (~ L*) galaxies at 2 ~ 5 —6 
(Capak et al. 2015, green pentagons) and 2 LBGs at 2 ~ 6 
(Willott et al. 2015, violet hexagons). Our model is in good 
agreement with the observations and local determinations 
(2 ~ 0) of the Lcii-SFR relation (De Looze et al. 2014), as 
we discuss in more details in Sec. 5. For simulated galaxies 
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brighter than Muv < —16, the Fpeak — Mw relation shows 
a bimodal trend: while 90% of the galaxies follow eq. 2, the 
remaining ~ 10% have fainter fluxes as a consequence of 
their lower metallicity valnes (see eq. 1). 

We finally note that we limit our comparison between 
simulation results and observational data to 2 ~ 6 nor¬ 
mal star forming galaxies (SFR 10^) without considering 
[C ii] observations in 2 ~ 6 QSO. In the latter case, obser¬ 
vations refer to galaxies hosting super massive black holes 
(M. ~ lO^M©) that are not present in the P14 simulation, 
because of its limited volume. 

2.3 CMB scattering from intergalactic metals 

The IGM and CGM are not seen in emission, since their typi¬ 
cal densities are too low (n < 0.1 cm“®) for the upper levels 
of the [G ii] transition to be efficiently populated through 
collisions with electrons and/or protons (Suginohara et al. 
1999; Gong et al. 2012). In this case the spin temperature 
of the transition approaches the GMB one'^. 

As a result a G ii ion at rest with respect to the CMB 
cannot show up in emission or absorption against this back¬ 
ground (e.g. de Bernardis et al. 1993; Maoli et al. 1996; da 
Cunha et al. 2013). However, if the ion has a peculiar ve¬ 
locity, for example approaching with a velocity Vp, it will 
receive a larger (i.e. Doppler-boosted) GMB photon flux in 
the direction of its motion. As CMB photons are resonantly 
scattered (e.g. Basu et al. 2004; Basu 2007; Schleicher et al. 

2008) by an atomic transition of frequency I'o — 1901 GHz 
the result is an emission feature in the GMB spectrum at 
i^obs = i^o/il + z). In the case of a receding C ii ion, the op¬ 
posite situation occurs, resulting into an absorption feature. 
G II ions in the IGM and GGM are therefore expected to 
generate GMB spectral distortions whose intensity depends 
on both their abundance and peculiar velocity, the latter 
being of the order of ~ 50kms“^. 

The differential amplitude A7j/ of the signal with re¬ 
spect to the CMB intensity is given by (Maoli et al. 1996): 

= (1 - (3 - a^)(vp/c) (3) 

where jAv)/ is the CMB spectral index and 

the optical depth (tj^) can be written as (e.g. Gallerani et al. 
2006; Pallottini et al. 2014b): 

= I'D (Tre^/m.er'o) NcWip , (4) 

where is the oscillator strength of the [C ii] transi¬ 
tion (e.g. Basu et al. 2004), e and me are the electron 
charge and mass, Acii is the [C ii] column density®, -i/i = 
'(p{(y — Vo)/A vd) is the normalized line profile (e.g. Meiksin 

2009) , Avd = {yolc)^j2kBTklmc is the thermal Doppler 

^ Gong et al. (2012) have shown that the IGM/CGM spin tem¬ 
perature of the [G ii] transition is close to the GMB temperature 
up to 2 ~ 2. At lower z, the soft UV background pumping effect 
may be sufficiently strong to decouple the spin temperature from 
the GMB temperature. Note that the redshift of the decoupling 
depend on the considered transition (e.g. Hernandez-Monteagudo 
et al. 2007, , for [O ij). 

® Note within our formalism the ISM does not contribute to Ti, 
in eq. 4, since ISM emission is accounted via the V15 model (see 
Sec. 2.2). 


broadening, ks the Boltzmann constant, me the carbon 
atom mass, and c the speed of light. We specify that we are 
focusing on CMB fluctuations produced by C ii ions, but 
the same formalism can be extended to other fine structure 
lines. 

In Fig. 2 we show an example of [C ii]-induced CMB 
distortions/fluctuations® for a single line of sight (l.o.s.) ex¬ 
tracted from a simulation snapshot centered at 2 = 6. In 
the left panel we plot Al^/Ii, as a function of frequency (v, 
lower axis) and velocity {v, upper axis); in the right panel, 
the metal patch responsible for the fluctuations is depicted 
by plotting the relevant physical characteristics (from top 
to bottom A, Z and Up) as a function of comoving distance 
(kpc/i“®) along the simulated l.o.s.. The origins of the ve¬ 
locity (left panel) and distance (right panel) axis are shifted 
for displaying purpose. 

Absorption {AIv/R < 0) and emission {AIv/Iv > 0) 
features are highlighted with a filled red and blue region, 
respectively. At i/ ~ 271.78 GHz we can see a single ab¬ 
sorption peak (labelled Ai in the Figure); the emission is 
instead characterized by a double peaked structure: Ei at 
V ~ 271.81 GHz and iJ 2 at i/ ~ 271.85 GHz. Both A\ and E\ 
have FWHMs of order ~ 20kms“®, while the FWHM of E 2 
is larger, ~ 40kms“^. The l.o.s. intersects the metal patch 
for ~ 0.2Mpc/i~®. The CMB interacts with a CGM peak 
(A ~ 10^) at a; ~ —30 kpc h~^ and an IGM peak (A ~ 10®'®) 
at a: ~ -1-75kpc/i~^. 

The GGM is responsible for both the absorption Ai 
and the emission E\ features. The change from absorp¬ 
tion to emission is caused by the change of sign in Vp. As 
Ncii oc ZA, the maximum of the GGM signal is produced 
rightwards of the A maximum, i.e. at a: ~ —5kpc/i~®. At 
this location Vp < 0, thus the GGM absorption feature (Ai) 
is stronger than the emission {Ep). The broader emission 
peak E 2 is caused by the IGM. Because of the lower density, 
this emission is one order of magnitude lower than the CGM 
absorption, and of the same order of the emission E\. 

The cumulative effect of the patch can be calculated by 
integrating over the shown bandwidth [Av = 1 GHz). This 
yields a net absorption, with intensity AI/I ~ —7 x 10“^. 
The order of magnitude of the cumulative effect is compara¬ 
ble with that of the amplitude of the fluctuations due to up- 
scattering of GMB photons with hot electrons (T > 10'* K) 
produced by first stars/galaxies radiation, supernova feed¬ 
back and structure formation shocks at redshifts 2 > 10 — 20 
(e.g. Chluba 2014). 


3 CMB FLUCTUATIONS MAPS 

We select a sample of ~ 4 x 10® l.o.s. extracted from a 
snapshot of the simulation at 2 = 6 to compute maps of 
CMB metal fluctuations. At this redshift, the P14 simulation 
is characterized by a maximum field of view (FOV) of ~ 
(350")^ and an adaptive spatial resolution down to ~ 500 pc. 


® In the CMB terminology, distortions indicate spectral varia¬ 
tions, while anisotropies refer to spatial variations. In this work 
we study spectral distortions that are induced by a non-uniform 
distribution of metals. Throughout the paper, we refer to them 
as fluctuations. 
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Figure 3. [C ll] column density map (left panel) and corresponding CMB fluctuations map (right panel) from a simulated field of view 
(FOV) of ~ (90”)^ at z = 6. Left panel: map of the [C ll] column density (A^cil) responsible for the CMB fluctuations. The column 
density is shown for the range 12 < log(AiQii/cm^) < 18 for displaying purposes. Right panel: the [C ll] CMB fluctuations map is shown 
as the brightness temperature AT = AIc^ / (2i'^kg) integrated on the simulation maximum available bandwidth, i.e. Au ~ 2.6 GHz. 
Note that the black spots inside the fluctuation patches (|AT| > 0) are caused by the ISM masking. In both panels, the angular scale is 
indicated as an inset. 


implying that maps are resolved at angular scales of 0res — 
0 . 2 ". 

Fig. 3 shows the [C ii] column density map (left panel) 
and the corresponding CMB fluctnation map (right panel) 
for a FOV ~ (90")^. The [C ii] colnmn densities are due to 
projection of the metal bubbles that originate around galax¬ 
ies, pollute the IGM/CGM and extend in the cosmic web fil¬ 
aments. The density profiles around galaxies are self-similar 
once scaled with the virial radius of the parent dark matter 
halo (Pallottini et al. 2014b), and metals usually extend out 
to ~ lOrvir. Thus, extended patches with high A^cii values 
are found nearby older (more massive) galaxies, which have 
more time to increase the metallicity of the central CGM 
part and pollute the surrounding IGM environment. As a 
reference, the ISM typically has log(Aicii/cm“^) 18, the 
CGM (IGM) can be enriched up to log(Aicii/cm“^) Si 16 
(log(Acii/cm“^) > 12). 

In the right panel of Fig. 3, CMB fluctuations are 
expressed in terms of differential brightness temperature, 
AT = Ale?/(2iPkB), where A7 is obtained by integrating 
All, (see eq. 3) on the total bandwidth of the simulation, 
i.e. Av ~ 2.6 GHz. The signal ranges from emissions up to 
AT ~ -l-lO^/rK down to absorptions of order AT ~ —lO^/rK. 
The signal is negligible (|AT| 10“^pK) in ~ 60% of the 

selected FOV, consistently with the analysis of the metal 
filling factor at this redshift (Pallottini et al. 2014a). 

The morphology of CMB fluctuations follows the Acii 
distribution. The lack of an exact match between the respec¬ 
tive maps is due to the dependence of AT from the peculiar 
velocity field. Absorption and emission features can arise 
from the differential velocity structure inside a single metal 


bubble. This is the spatial analog of the emission/absorption 
features as a function of v for the l.o.s. shown in Fig. 2. 

Fluctuations are preferentially located in correspon¬ 
dence of log(Acii/cm“^) > 12 patches, characterized by 
typical sizes 9 ~ 20"; local maxima of the fluctuations 
(|AT| > 10^/iK) are usually found in smaller spots {9 ~ 1") 
characterized by log(Acii/cm“^) ~ 17. This result is con¬ 
firmed through the analysis of the CMB fluctuations power 
spectrum that peaks at angular scale 0 ~ 1", as detailed in 
App. A. 

CMB fluctuations on large (~ 10") and small (~ 1") 
scales arise from the enriched IGM and GGM, respectively. 
Smaller scales correspond to the ISM of galaxies. These re¬ 
gions have been masked^, since they do not contribute to 
CMB fluctuations, while possibly shining as [C ii] emitters. 

The [C ii] is a biased tracer of the cosmic web as it 
is mostly concentrated in the intersection of the filaments 
(“knots” hosting galactic environments) where the density 
is sufficiently large (A Si 10^) that recombinations overcome 
photoionization from the UV background (and possibly lo¬ 
cal) radiation field. In addition, as several observational and 
theoretical works have shown, the filling factor of metals at 
these high -2 is relatively small, of the order of 10%-20%. So 
independently of the ionization state, it would be anyway 
difficult to find carbon in regions of moderate to low density 
(A 5) corresponding to the filaments. 


^ The ISM masking interests less than 0.5% of the total FOV. 
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4 MOCK OBSERVATIONS 

Starting from our model for [C ii]-induced CMB fluctuations 
and [C ii] emission from z = 6 galaxies, we generate mock 
observations specifically suited for comparison with ALMA 
data. The [C ii] line at 157.74 /rm from 2 = 6 is redshifted 
in ALMA BAND6 (1.1 — 1.4mm). At this wavelengths the 
FOV (primary beam) of ALMA antennas is ~ 20”, and the 
angular resolution (synthesized beam) ranges from 1.8” to 
0.14”, for the most compact and the most extended conhg- 
uration, respectively (e.g. Maiolino 2008)®. 

We account for ALMA angular resolution by convolving 
our maps with the synthesized beam, that can be character¬ 
ized as a 2D-Gaussian. In this Sec., we assume a major axis 
of 1 . 4 ”, a minor axis of 0.9”, and a position angle of 91.6 deg. 

We note that the synthesized beam that optimizes the 
intensity of the signal depends on the experiment that one 
would like to carry on. In the case of CMB fluctuations, the 
synthesized beam that optimizes the chance of detection is 
~ 1 ”, since the power spectrum peaks at this scale, as shown 
in Sec. 3 and by the power spectrum analysis, reported in the 
App A. Viceversa, for detecting [C ii] emission from galaxies, 
a synthesized beam of ~ 0.5” is more appropriate since it 
corresponds to typical galaxy sizes at 2 ~ 6, thus minimizing 
the surface-brightness bias for galaxy detection (e.g. Grazian 
et al. 2011). This different observational set-up is further 
described in Sec. 5. 

GMB fluctuations and galaxy emission maps are gener¬ 
ated by integrating on Azz = 8 GHz, i.e. the maximum band¬ 
width available for continuum observations. This is done 
in order to maximize both the sensitivity and the possi¬ 
bility of finding GMB signal in the selected FOV. The P14 
simulation size limits the maximum available bandwidth to 
Au = 2.6 GHz (see Sec. 3). We generate a longer light-cone 
by replicating the 2 = 6 simulation box. To avoid spuri¬ 
ous periodicity effects, we randomize the simulation box by 
applying translations, rotations and reflections (e.g. Blaizot 
et al. 2005). Then, the signal from the light-cone is inte¬ 
grated over Au = 8 GHz, by taking into account that ALMA 
continuum band is not contiguous, i.e. it is composed of two 
4 GHz bands, separated by 10 GHz. 

Finally, mock maps are combined with noise maps gen¬ 
erated with ALMAOST (Heywood et al. 2011). Within the pri¬ 
mary beam, the noise map has a Gaussian distribution with 
zero mean and r.m.s. crjv; the r.m.s. accounts for instrumen¬ 
tal uncertainties and atmospheric conditions. We assume 
0.913 mm of water vapor for our fiducial observational set¬ 
ting. Using the ALMA full array, we have a continuum sensi¬ 
tivity of CTjv — 11.09 (t/hr)~^/^(Azz /8 GHz)~^'^^[^Jy/beam] 
for an observational time t at Vahs — 271 GHz. 

In Fig. 4 we plot mock maps centered on the same 
Muv — —19 galaxy. In the upper panels, we plot separately 
the galaxy emission (left) and CMB fluctuations (right) sig¬ 
nals, expressed in /rJy/beam. The angular scale and the 
synthesized beam are indicated in the inset. In the galaxy 
emission map, the symbol ‘X’ indicates locations of galax¬ 
ies, which are labelled according to their Muv- Galaxies are 
mostly found clustered in the upper region, and the plot 

® For further details, please refer to the ALMA handbook 
https://almascience.eso.org/documents-and-tools/cycle-O/ 
alma-technical-handbook/at_download/file 


shows a clear correlation between galaxy UV magnitude and 
[C ii] emission (see the right panel of Fig. 1 in Sec. 2.2). In 
fact, brighter [C ii] emitters {Al ~ /rjy/beam) correspond 
to brighter UV galaxies (Muv — —19), characterized by 
a total [C ii] flux of the order of Al ~ 100 /rJy; dimmer 
galaxies (Muv — —16.5) are identihed through [C ii] emis¬ 
sion peaks of the order of A7 ~ 10“^/rJy/beam; the faintest 
galaxies in the map (Muv — —15) show very low [C 11 ] emis¬ 
sion, Al 10“^^Jy/beam. 

The CMB fluctuations map (upper left panel of Fig. 4) 
shows various peaks both in absorption and emission. The 
signal due to CMB fluctuations result to be more extended 
and offset with respect to galaxy emission (indicated as ‘X’). 
This is expected from our analysis of Vcii and CMB fluctu¬ 
ations maps (see Fig. 3). Typically, emission/absorption fea¬ 
tures range between 10“^ |A/|/(/rJy/beam) 10“^, with 

local intensity maxima comparable to the emission signal of 
a Muv — —17 galaxy. The strongest peak occurs in absorp¬ 
tion and it characterized by A7 ^ — 10“^/rJy/beam. 

In the lower panels of Fig. 4, we plot the total signal, i.e. 
including the contribution from CMB fluctuations, galactic 
emission and two different noise levels: crjv — 0.1/rJy/beam 
(left) and un — 0.02/iJy/beam (right) that correspond to 
t ~ 1.3 X lO'^hr and t ~ 3.1 x 10®hr observing time, respec¬ 
tively. In these panels, the intensity of the signal is rescaled 
to the noise r.m.s.; additionally, detections of CMB fluctu¬ 
ations and galaxy emission are indicated with C; and Gi, 
respectively. 

For (j]v — 0.11/iJy/beam, the noise covers most of the 
CMB fluctuations and all the faint peaks of the galactic 
emission map {Muv ~ — 17). Detectable signals arise from 
the two Muv — —19 galaxies (Gi and G 2 ) and from the peak 
of the CMB fluctuations (Ci); the galaxies show > lOuiv 
emission peaks, and the fluctuation appears as a ~ 7ctjv 
absorption signal. 

In the case of the deeper mock observation, for (Tjv — 
0.02 /rjy/beam, Muv — —17.6 galaxies are detectable (G 3 
and G 4 ) at > bcrv as [C ii] emitters within the primary 
beam. Note that G 3 is a > lOuiv peak, while G 4 show an 
emission ~ 6 (Jjv, as its signal is partially suppressed by the 
nearby absorption from CMB fluctuations at Gi. Three more 
CMB fluctuations peaks are visible on the map: a ~ huiv 
absorption arising outside the primary beam (G 2 ), and two 
extended ~ 5 (tjv emissions (G 3 and G 4 ) located within 2.5” 
from the Muv — —19 galaxies. 

Note that CMB fluctuations in emission can be mis¬ 
taken for signals arising from galaxies. In principle, the 
origin of the signal can be distinguished by analyzing its 
spectral properties. Typically, galaxy [C ii] emission lines 
have FWHMs ^5 50 — 100kms“^ (e.g. Vallini et al. 2013; 
Decarli et al. 2014; Cormier et al. 2015), while CMB fluc¬ 
tuations usually show a double-peak structure in the spec¬ 
tra (see Sec. 2.3, in particular Fig. 2), with narrower lines 
(FWHM~ 30 —40kms“^) both in absorption and emission. 

A more robust discriminant for the nature of the signal 
can be provided by multiwavelength analysis. By compar¬ 
ing the [C ii] map with UV imaging of the held, we can dis¬ 
tinguish between the emission from CMB huctuations and 
galaxies. Knowing the location and UV luminosity of galax¬ 
ies in the held, we can infer the expected [C ii] hux (see Sec. 
2.2, in particular right panel of Fig. 1). JWST will have a 
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Figure 4. Mock ALMA [C ii] continuum {Au = 8 GHz) map observed with BAND6 with a synthesized beam of ~ The maps 

are shown for a total FOV of (38.04 ~ ns)^, twice BAND6 primary beam. In each panel, the map, the angular scale and the synthesized 
beam are indicated as insets. In the upper panels we separately plot the contribution for galaxy emission and CMB fluctuations. In the 
lower panels two mock observations (galaxy emission and CMB fluctuations and noise) for two different noise levels are shown. Upper 
left panel: [C ii] signal from galaxy emission. The positions of galaxies are indicated with an ‘X’ and labelled with the corresponding 
UV magnitude (Afuv? details in the text). The intensity is expressed as //Jy/beam and, to better appreciate the presence/absence of the 
signal, the intensity is in symmetric logarithmic scale. Upper right panel: [C ii] signal from CMB fluctuations from the same field of 
view. As for emission from galactic ISM, the intensity scale is nonlinear. Lower left panel: mock observation, obtained by summing the 
contributions from CMB fluctuation, galaxy emission and noise. The noise distributions fit a Gaussian with r.m.s. — 0.11/iJy/beam. 
To augment the readability of the map, the intensity is rescaled with (tn and the scale is limited to |A/| < lOcrjv- To help the description 
in the text, detections of CMB fluctuations and galaxy emission are indicated with Ci and Gi, respectively. For reference, we plot 
the ALMA primary beam and twice the primary beam with dashed lines. Lower right panel: as in the lower left panel but with 
(Jtv ~ 0.02/iJy/beam. 
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reference limiting magnitude of Muv — —16 at 2 = 6 , thus 
it seems to be perfectly suited for such comparison. 


5 OBSERVATIONAL PROSPECTS 

The amplitude of [C ii]-induced CMB fluctuations is ex¬ 
pected to be |A/| ^ 1/iJy/beam, as discussed in the pre¬ 
vious section. With the aim of defining the optimal obser¬ 
vational strategy for detecting such weak signals, we first 
consider the feasibility of a blind survey experiment with 
ALMA. We compute the probability distribution functions 
(PDFs) of the intensity of [C ii]-induced CMB fluctuations 
and [C ii] galaxy emission at 2 = 6 . We consider the sig¬ 
nals extracted from the whole simulation FOV, ~ (350”)^, 
on the Ai/ ~ 2.6 GHz bandwidth centered at the frequency 
z^obs — 271 GHz. For both signals, the intensity accounts for 
the ~ (1”)^ ALMA synthesized beam introduced in Sec. 4. 
For each PDF, the corresponding integral is normalized to 
unity. 

In the left panel of Fig. 5, we plot the PDFs for 
CMB fluctuations and galaxy emission with blue solid and 
red dashed lines, respectively. The plot is cut at |A/| = 
10“^/iJy/beam for displaying purpose. As a reference, with 
vertical lines we plot 1-a flux detection limits that corre¬ 
spond to the ALMA observing time indicated in the leg¬ 
end, calculated with the ALMA Sensitivity Calculator®. 
For the CMB fluctuations, the PDF rapidly decreases from 
10“^'® at |A/| ~ 10“®^Jy/beam down to 10“'* at |A/| ~ 
10“® ®/rJy/beam. For the galaxy emission, the PDF de¬ 
crease is less steep, and the two PDFs cross at |A7| ~ 
10“*/iJy/beam. 

We indicate with ctn the 1-a flux detection limit of 
the ALMA blind survey. Then the integral of the PDF 
above ctn provides the probability Pdis and Pgai of detect¬ 
ing CMB fluctuations and galaxy emission, respectively. For 
Ctn = 10“®^Jy/beam (corresponding to an ALMA observing 
time t ~ 10®hr), we find Pdis — 4% and Pgai — 1%. Given the 
current telescope sensitivity, if we consider a more realistic - 
thought challenging - observation, i.e. ctn = 10“*/iJy/beam 
{t ~ 10 *hr), we find Pgai cs Pdia — 0.05%. This highlights 
that it is extremely difficult to detect GMB fluctuations with 
a blind survey. 

Therefore, as an alternative strategy we consider the 
search for CMB fluctuations in the close proximity of high- 
2 galaxies. We compute the cross-correlation R) between 
CMB fluctuations peaks and UV galaxy emission. We cal¬ 
culate 5 by using an extension of the Hamilton estimator 
(e.g. Nollenberg & Williams 2005): 

aO) + 1 = {DGDc{DRDn))/{{DGDR){DcDR )), (5) 

where DiDj is the number of data pairs from sets i and 
j within angular distance 9, the subscripts G and C indi¬ 
cate the set of positions for galaxy and CMB fluctuations 
peaks respectively, and R labels a set of locations extracted 
from a random uniform distribution; the operator (...) in¬ 
dicates the average on different realizations of R. Both the 
dimension of the set R and the number of realizations for 
the average are fixed to reach a suitable convergence for 

® https://almascience.eso.org/proposing/ 
sensitivity-calculator 


For what concerns CMB fluctuations, we consider only 
strong peaks, i.e. fluctuations with |A/| > O.l/iJy/beam. 
Within the simulated FOV we find ~ 30 of such peaks, with 
signal intensity up to 0.25/iJy/beam in emission and down to 
—0.82/iJy/beam in absorption. For these signals, the mean 
|A/| is ~ 0.23/iJy/beam and the r.m.s. is ~ 0.25/iJy/beam. 
For what concerns galaxies, we select sets of objects charac¬ 
terized by different UV limiting magnitudes. 

In the right panel of Fig. 5, we plot with red, blue and 
black lines the cross-correlation calculated for galaxies with 
Afuv < —18, —16 and—14, respectively. Although there is a 
clear cross-correlation between strong CMB fluctuations and 
galaxy emission, we find no strong dependence of the cor¬ 
relation scale from the considered Muv limits*®. Therefore, 
the most promising - though challenging - observational 
strategy for detecting [C ii]-induced CMB fluctuations is to 
perform a deep (ctjv O.l/rJy/beam) ALMA BAND6 obser¬ 
vation pointing at a Muv — —19 known galaxy, e.g. in the 
HDF south, where the deepest photometry is available. This 
observational strategy is detailed in Sec. 4. 

From Fig. 5 it appears that [C ii] emission detection is 
very challenging even from the bright end galaxies in our 
simulation (Muv ~ — 19). As anticipated in Sec. 4, this 
is critically dependent on the assumed ~ 1” synthesized 
beam. A galaxy at Muv — —19 has a [C ii] emission of 
~ lOO^Jy (see Fig. 1 and eq. 2). This correspond to a flux 
of ~ 1/rJy/beam in the ALMA map with a ~ 1” synthesized 
beam. Hence, with this set-up, ~ 10® hours of integration 
time are required in order to detect the signal at ~ 4 ct. 

The detectability drastically improves when using in¬ 
stead a synthesized beam similar to the galaxy size, i.e. 
~ 0 . 5 ”, like the one used in recent high -2 observations 
(Maiolino et al. 2015; Capak et al. 2015; Willott et al. 2015). 
With this set-up, by using eq. 2 and the ALMA Sensitivity 
Calculator, we find that a Muv = —18, —19, —20 galaxy can 
be detected at 4 ct in ~ 2000, 40,1 hours, respectively. Indeed, 
[C ii] has been detected in high -2 galaxies with magnitudes 
up to Muv > —20 with f ~ 1 hr (Maiolino et al. 2015; 
Willott et al. 2015). With a more demanding observational 
time of ~ 40 hr, we would be able to detect Muv — —19 
galaxies at 2 ~ 6, thus sampling galaxies more similar to 
the true reionization sources. 

As noted in Pallottini et al. (2015), at 2 = 6 P14 galax¬ 
ies have a specific SFR of sSFR = SFR/M^ ~ 2.5Gyr“*, 
that is consistent with observations by Daddi et al. (2007) 
and Gonzalez et al. (2011). However, more recent observa¬ 
tions seem to favor an higher sSFR (e.g. Stark et al. 2013, 
sSFR = 5 Gyr“*). The latter sSFR would yield a [C ii] emis¬ 
sion higher by a factor ~ 2, namely a required observing 
time a factor ~ 0.7 smaller. An increase of sSFR results in 
an increase of metal production, that may enhance the sig¬ 
nal of CMB fluctuations by up to a factor ~ 2. However, 
it is difficult to quantify this enhancement, because of the 
complexity of the metal enrichment process. 

*® The cross-correlation scale does not strongly depend on the 
considered Muv limits because of the following. Brighter (more 
massive) galaxies are expected to be surrounded by more ex¬ 
tended and enriched CGM. However CMB fluctuations linearly 
depend on the velocity field, thus there is a lack of correspon¬ 
dence between high Nqh values and strong fluctuations (see Sec 
3 for details, in particular Fig. 3). 
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Figure 5. Left panel: PDFs of the absolute value of the signal from CMB fluctuations (blue solid line) and galaxy emission (red solid) 
are plotted. Both signals account for the assumed ALMA synthesized beam (see Sec. 4). The PDF integral is normaiized to unity, and 
the plot is cut at |A7| = 10“^/rjy/beam for displaying purpose. As a reference, with vertical lines we plot several flux detection limits: 
the corresponding ALMA observational time is indicated in the iegend. Right panel: Angular cross-correlation (^) function between 
high CMB fluctuation peaks (|A/| > 0.1/rJy/beam, see text) and gaiaxies with limiting magnitude Mpv = ~18i ~16 and ~ 14 (plotted 
with red, biue and biack lines, respectively). 


6 CONCLUSIONS 

We have studied the possibility of mapping heavy ele¬ 
ments via far infrared (FIR) emission from the interstel¬ 
lar medium (ISM) of high-z galaxies and cosmic microwave 
background (CMB) fluctuations induced by metals in the 
intergalactic medium (IGM). We focus on high -2 low mass 
(M* IO^^Mq) galaxies, which are expected to be abun¬ 
dant at 2 ~ 6 , and to represent the first (2 ^ 10 ) efficient 
metal polluters of the IGM (e.g. Madau et al. 2001; Fer¬ 
rara 2008) and sources of reionization (e.g. Barkana & Loeb 
2001 ). 

Among the FIR emission lines, the 
[G II] (^^3/2 ^ ^Pi ^ 2 ) transition at 157.74 is the 

brightest (e.g. Crawford et al. 1985; Madden et al. 1997), 
and it is used to trace and characterize 2 ~ 6 galaxies (e.g. 
Carilli & Walter 2013). We compute both the [C ii] emission 
arising from the ISM of 2 = 6 galaxies and the amplitude of 
[C ii]-induced^^ CMB fluctuations from metals in the IGM 
and circumgalactic medium (GGM). 

We use state-of-the-art high -2 hydrodynamical simula¬ 
tions (Pallottini et al. 2014a) that follow the evolution of 
the cosmic metal enrichment from 2 = 10 from 2 = 4. We 
calculate the galactic [G ii] signal by using an ISM model 
(Vallini et al. 2013, 2015) that accounts for the detailed sub- 
kpc scales structure of the emission, namely for the cold and 
warm neutral medium and in high density photodissociation 
regions (PDRs). We calculate the CMB fluctuations induced 
by C II ions by modelling the resonant scattering between 
the CMB and metals in the IGM/CGM (e.g. Maoli et al. 
1996; Basu et al. 2004) and we predict the related emission 
and absorption features arising in the FIR band. 

We then carefully analyze the theoretical signal of CMB 

The presented method can be straightforwardly extended to 
other FIR lines, as [N ll] at 122 pm, [O III] at 88 pm, [Si l] at 
129 pm and [O l] at 63 pm. 


fluctuations from the IGM/CGM. While the effects of res¬ 
onant scattering of GMB photons by very early metals 
(2 ^ 10 ) were studied with analytical models and simple 
assumptions about the metal distribution (e.g. Basu et al. 
2004; Schleicher et al. 2008), the present study is the first 
based on hydrodynamical numerical simulations extending 
to observable epochs (2 ~ 6 ). In terms of the differential 
brightness temperature, AT = A7c^/(2iz^fcs), we calculate 
the signal from 2 = 6 integrating on the simulation available 
bandwidth, Aiz = 2.6 GHz. We find that the metal-induced 
fluctuation signal is in the range AT ~ ilO^pK, i.e. can be 
seen either in emission or absorption. The peak of the sig¬ 
nal is found on scales of 0 ~ 1”, in correspondence of CGM 
absorption systems characterized by C ii column density of 
log(Acii/cm“^) ~ 16. 

To test the detectability of the [C 11 ] signal, we have con¬ 
structed and analyzed mock observations specifically suited 
for comparison with ALMA BAND 6 data. 

We predict that [C ii] emission is correlated with Muv- 
At Muv < — 20 , the faintest high -2 galaxy from which 
[G ii] emission is detected, our relation (eq. 2) is in good 
agreement with recent observations (Maiolino et al. 2015; 
Capak et al. 2015; Willott et al. 2015). We find that a 
Muv = —18, —19, —20 galaxy can be detected at 4cr in 
~ 2000, 40, 1 hours, respectively. Indeed, [C ii] has been de¬ 
tected in high -2 galaxies with magnitudes up to Muv > —20 
with t ~ 1 hr (Maiolino et al. 2015; Willott et al. 2015). 
With a more demanding observational time of ~ 40 hr at 
2 ~ 6 we would be able to detect Muv — —19 galaxies, 
thus sampling galaxies more similar to the true reionization 
sources. Additionally, the relation is the analogous of the 
local (2 = 0) Lcii-SFR relation (De Looze et al. 2014). Eq. 
2 highlights the possibility of using [C ii] as a tracer of star 
formation activity, that can be used independently from UV 
flux determination. This is particularly important at high -2 
as, contrary to UV light, [C 11 ] emission is not affected by 
dust extinction. 
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FIR observations allow in principle to simultaneously 
detect galajdes and to map their surrounding CGM through 
[C ii] induced CMB fluctuations. The detection of this sig¬ 
nal would be a breakthrough for our understanding of the 
early phases of galaxy formation and cosmic enrichment. 
The CGM is the interface regulating the outflows of enriched 
material from the galaxy and the inflows of pristine gas from 
the IGM; it has been probed so far up to 2 ~ 2 with absorp¬ 
tion line experiments towards background sources, QSO or 
galaxies (e.g. Steidel et al. 2010; Churchill et al. 2013; Liang 
& Chen 2014). The intervening CGM associated with a fore¬ 
ground galaxy leaves an absorption feature in observed spec¬ 
tra. With a large number of galaxy-absorber pairs, it is then 
possible to statistically determine the equivalent width of a 
given absorption line as a function of the line of sight impact 
parameter. Detecting CMB fluctuations is equivalent to ac¬ 
tually map the bidimensional distribution of metals in the 
CGM, and it would allow to extend the enrichment study 
at higher redshift. 

We have compared the efficiency of different observa¬ 
tional strategies for detecting [C ii] induced CMB fluctu¬ 
ations. We find that this signal is very faint, e.g. |A/| ~ 
10“^/iJy/beam. Because of the signal amplitude and the 
[C ii] emission filling factor, we find that a blind experi¬ 
ment yields a low probability for the signal detection. How¬ 
ever, in our maps we see strong (i.e. j A/j 10“^/rJy/beam) 
CMB fluctuations, that are characterized by a mean |A7| 
of ~ 0.23/iJy/beam a r.m.s. of ~ 0.25/rJy/beam and maxi¬ 
mum |A7| up to 0.82/iJy/beam. On average, these fluctua¬ 
tions have signals comparable to the typical [C 11 ] emission 
from a Muv — —18 galaxy. Our analysis highlights that 
strong CMB fluctuations are found typically within ~ 10” 
of galaxies, regardless of their Muv- 

We have also considered an alternative observational 
strategy to detect (strong) CMB fluctuations that consists of 
a deep ALMA pointing in a field where known Muv — —19 
galaxies (e.g. in the HDF south) are present. Our model pre¬ 
dicts that, for an ALMA sensitivity of av = O.l/rJy/beam 
CMB fluctuations would be detected at a confidence level 
a > 3, 7 with a probability ~ 25%, 5%, respectively, 
while for ajv — 0.05/rJy/beam we expect a detection of 
CMB fluctuations with c.l. cr > 3, 5, 7 with a probability 
~ 70%, 30%, 15%, respectively. We note that the sensitiv¬ 
ity requested for detecting [C 11 ] induced CMB fluctuations 
(cTjv ,$ 0.1/rJy/beam) corresponds to extremely long ALMA 
observing times (t ^ lO^hr), making the detection of this 
signal extremely challenging with current facilities. 

However, by stacking deep ALMA observations (i/obs — 
272 GHz) of several Muv — —19 (lensed) galaxies it may be 
possible to “statistically” detect such elusive CMB fluctua¬ 
tions induced by C ii ions at 2 ~ 6. The same experiment, 
repeated at different frequencies (z^obs — 1900/(1 -I- 2 ) GHz) 
would allow to track the metal enrichment history through 
cosmic times. We defer this further statistical analysis to a 
future study. 
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APPENDIX A: POWER SPECTRUM 

We investigate the CMB fluctuations morphology by calculat¬ 
ing the angular power spectrum (PS, Ci) of their amplitude. At 
2 = 6 the simulation FOV is limited to ~ (350'')^, therefore 
the minimum angular scale solved is Z ~ 10^. This allows us to 
compute the Ci in the fiat-sky approximation (e.g. White et ah 
1999; Mesinger et ah 2012). We calculate PS for different band- 
widths 20 MHz ^ Az/ ^ 2.6 GHz, where the upper limit of the 
interval represents the total bandwidth of the simulation. For 
;S 1.3 GHz, this procedure provides multiple bandwidths that 
we assume to be independent to calculate the averaged Ci . 

In upper panel of Fig. Al, we plot the power spectrum as 
Pi = l(l-\-l)Ci/2n as a function of the angular scale (9 = 27r/Z) for 
different bandwidths, as indicated in the legend. Independently 
of Ah', the power spectrum has a noise-like shape, i.e. Pi oc Z^, 
for Z ^ 10®. In correspondence of smaller scales. Pi reaches a peak 
at 0 ~ 1 " and then flattens, confirming the qualitative analysis 
discussed in Sec. 3. 

In the central panel of Fig. Al, we compare our result (Ai^ ~ 
2.6, solid black line) and PS for CMB distortions/fluctuations in¬ 
duced by primordial chemistry (Basu 2007, violet dotted line). 
The latter PS is obtained for C ii that are assumed to be uni¬ 
formly distributed (A = 1) at^^ 2 = 4 with Z = 10~^Zq,. With a 
grey shaded region we indicate the cosmic variance of the [C ii]- 
induced CMB fluctuations. As a reference, we plot the primary 
CMB power spectrum inferred from Planck observations (Planck 
Collaboration et al. 2013, red crosses) and calculated with CMB- 
FAST (Zaldarriaga & Seljak 2000, green dashed line). 

For the primordial fluctuations. Pi has a functional depen¬ 
dence on Z as the primary CMB power spectrum (Basu 2007; 
Schleicher et al. 2008). With respect to the primary power spec¬ 
trum, these induced fluctuations produce a ~ 10^ smaller power, 
e.g. Pi ^ fiK^ on scales Z ^ 10^ (Basu et al. 2004; Basu 2007). 

We define, the PS peak as (AT) = ^/Pi at 0 ~ 1" ^ and we 
study its dependence on Av. The result is shown in the lower 
panel of Fig. Al, where we plot the PS peak mean value (red 
circles), the r.m.s. variance (green pentagons) and the max/min 
values (blue triangles) as a function of Ai^. We restate that the 
mean and variance are calculated for the multiple bandwidths 
extracted from the simulation: they must not be confused with 
variation of CMB fluctuations from different metal bubbles. 

The peak value increases with decreasing bandwidth up to 
Au ~ 300 MHz, and decreases for smaller Ai^. This behavior can 
be explained as follows. The peak value increases as the band¬ 
width becomes comparable to metal bubble size in frequency 


As shown in Basu (2007, in particular, see left panel of Fig. 
1), at a different redshift, the primordial chemistry PS induced 
by C II is comparable to the one reported here. 
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space (see Fig. 2). As we further decrease Az/, it becomes increas¬ 
ingly difficult to find enriched structure in the selected bandwidth, 
and the peak becomes shallower. 

Finally, note that the peak (y/Pi) increases faster than 
for 2.6 ^ Ai//GHz ^ 0.3 and decreases for smaller Az/. The power 
spectrum of a pure noise is expected to behave as Az/~^/^. There¬ 
fore - in principle - the peculiar trend with Az/ can discriminate 
a signal from CMB fiuctuations from noise. 



Figure Al. [C ii]-induced CMB power spectrum (PS, Pi = 
I {I 1) Cil27r) as a function of the angular scale, I and 6 in the 
lower and upper axis, respectively. In the upper panel we plot 
Pi calculated for different bandwidths, as indicated in the legend. 
In the central panel we compare our result with other PSs, as 
detailed in the legend and in the text. A color version of the 
plot is available in the online version of the paper. Lower panel: 
maximum of PS ((AT) = y/Pi at 0 ~ 1") as a function of the 
bandwidth (Az.^). For each Ai^ we plot the mean (red circles), the 
r.m.s. variance (green pentagons) and the max/min values (blue 
triangles). The maximum is indicated in unit of fiK (left axis) 
and /ijy (right axis). 
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